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ABSTRACT 

The time-dependent Comptonized output of a cool soft X-ray source drifting 
inward through an inhomogeneous hot inner disk or corona is numerically simulated. 
We propose that this scenario can explain from first principles the observed trends in 
the hard time lags and power spectra of the rapid aperiodic variability of the X-ray 
emission of Galactic black- hole candidates. 



Subject headings: X-rays: stars — accretion, accretion disks — black hole physics - 
radiative transfer — radiation mechanisms: thermal 



Introduction 



The rapid aperiodic variability of the X-ray emission of Galactic black-hole candidates 
(GBHCs) and low-mass X-ray binaries (for a review see van der Klis 1995| ) is attracting an 



increasing amount of attention since spectral information alone is generally not sufficient to 
distinguish between different models for the high-energy emission of X-ray binaries. Its power 
spectrum density (PSD), autocorrelation function, time lags between different energy bands, 
coherence etc. provide important diagnostics of important parameters such as the extent of the 
emitting region, the dominant microscopic timescales, the dominant emission mechanism and the 



geometry of the source (Kazanas et al. 1997, Hua et al. 1997, Bottcher & Liang 1998) 



Early measurements of the Fourier-frequency dependence of time lags between the signals in 
different X-ray energy channels from GBHCs ( Miyamoto et al. 1988] , 1993) have been interpreted 



as evidence against Comptonization of soft photons in a hot, uniform plasma as the mechanism for 
the production of the rapidly variable hard X-rays. Recently, Kazanas et al. (1997) and Hua et al. 
(1997) have pointed out that Comptonization of internal soft photons in an inhomogeneous hot 
corona can reproduce the observed PSD slopes and hard time lags of Cyg X-l and other Galactic 
black-hole candidates in the low-hard state. Bottcher & Liang (1998) have expanded this idea 
further and discussed observable signatures to distinguish between different source geometries. 
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However, Bottcher &; Liang (1998) also pointed out that any scenario in which the observed 
hard time lags are purely due to static Comptonization, requires that the radial extent of the 
hot corona exceeds ~ 10 4 Schwarzschild radii of a solar-mass black hole (R ;> 10 10 cm). This is 
incompatible with current models of accretion flows onto Galactic black holes (for a recent review 
see Liang 1998| ) even from simple energy arguments. Moreover, according to Vaughan & Nowak 



(1997) a pure Comptonization scenario might have problems reproducing the observed coherence 
functions due to the stochastic nature of the Comptonization process, and Hua et al. (1997) 
mention that the coherence observed in Cyg X-l can only be reproduced by Comptonization in an 
inhomogeneous hot corona if the physical conditions of the Comptonizing region remain constant 
over long time scales (hours) . This suggests that additional variability of the physical parameters 
of the Comptonizing region may be involved in the rapid aperiodic variability of GBHCs. 

A first model along these lines was developed by Poutanen &: Fabian (1998) who suggested 
that the Comptonization occurs in active regions in an optically thin corona above a cold accretion 
disk ( [Haardt fc Maraschi 199l| , 1993, (Stern et al.T995|) . The active regions were assumed to be 



heated, e. g., due to magnetic flares ( Nayakshin fc Melia 1997| ), causing the coronal temperature 



to rise exponentially with time before it drops off sharply. With these temperature profiles, the 
observed X-ray spectra and hard time lags of Cyg X-l could be reproduced. However, the flare 
time profile, in particular the instantaneous temperature drop at the end, required to yield the 
observed hard time lags, seems rather artificial. There also seems to be a problem with the 
coherence function in such a scenario since it would require that the entire corona flares uniformly 
and homogeneously. 

In this Letter, we propose an alternative and, we believe, more natural model in which 
a blob of cool, dense matter spirals inward through an inner hot corona until it disappears 
at the event horizon of the black hole. The inner portions of the accretion flow are assumed 
to form a quasi-spherical, inhomogeneous, hot, optically thin corona, which could either be a 
Shapiro-Lightman-Eardley (1976) type accretion disk or an advection dominated or transonic 
accretion flow ( [Narayan fc Yi 1994| , |Chen et al. 1995| , [Esin et al. 1997| , [Luo & Liang 1998j ). This 



picture is partially motivated by the 2-phase disk scenario recently proposed by Krolik (1998). 
As the cool, dense blob drifts inward, the soft X-ray radiation it emits is Comptonized in an 
increasingly hot and dense corona, which results in the Comptonized X-ray spectrum becoming 
harder with time, causing the observed hard time lags. The key idea here is that the observed 



cutoff period of seconds in the hard-time-lags versus Fourier-period plots (e. g., Miyamoto et al. 
1988, Cui et al. 1997a| ) is associated with the radial drift time scale of the blob. 



Details of the numerical scheme adopted to simulate the above scenario as well as an extensive 
parameter study will be published elsewhere (Bottcher & Liang 1999, in preparation). The 
Monte-Carlo code used here has been described in Bottcher & Liang (1998). In Section 2, we 
describe the model and outline the method of our numerical simulations. Resulting power spectra 
and time lags are presented and discussed in the light of observational results in Section 3. Section 
4 contains a short summary and our conclusions. 



- 3 - 



2. The model 

The basic model presented here consists of a soft X-ray source (the "blob"), emitting a 
thermal blackbody spectrum of temperature Tbb with 0.1 keV ^ kTsB ^ 1 keV, which is drifting 
inward with a constant radial drift velocity /3 r c. Tbb is allowed to have a power-law dependence 
Tbb oc r~ PBB on the distance r from the center of the black hole. When the soft photon source 
reaches the radius rj n of the last stable orbit of the accretion disk, it rapidly falls onto the event 
horizon of the black hole and disappears. In the inner ~ 100 Rs around the black hole, we assume 
the existence of a quasi-spherical, hot, inhomogeneous corona. For the simulations presented in 
this Letter, we fix the radial dependence of the density n c and temperature T c in this central 



corona to values representative of an advection dominated accretion flow (Narayan & Yi 1994 



Chen et al. 1995| ), i. e., a temperature profile T c oc r 1 and a density profile n c oc r 3 / 2 are 



assumed. 

To calculate the observable spectra and temporal variability resulting from this model 
we perform time-dependent Monte-Carlo simulations of Comptonization in an inhomogeneous 
medium. The Comptonization Code we use is based on the code developed by Canfield et al. 
(1987, see also [Liang 1993| and Bottcher fc Liang 1998 ). The radial profile of the corona is 



approximated by dividing the corona into ~ 50 radial zones of approximately equal Thomson 
depth increments. Within each zone the density and temperature have a constant value. The 
inner boundary of the corona is treated as a purely absorbing boundary (the black hole horizon). 
We are simulating a single blob drifting inward through this corona. The photons escaping the 
corona at its outer boundary r ou t ~ 100 Rs are sampled in small time bins (At <^ 0.01 s) and in 5 
energy bins. The resulting light curves in all energy bins are then Fourier transformed, and the 
corresponding power spectra and hard time lags are calculated. 

A rough analytical estimate of the mean energy of escaping photons as a function of escape 
time can be found based on the mean number of scatterings, N sc (r em ) ~ maxjry, Tj.}, where tt 
is the radial Thomson depth of the corona from the point of emission of a soft photon at radius 
r em , to the outer boundary. Due to the steep density gradient, most of the Comptonization 
will take place around r em . Thus, the escaping photon energy is E esc m E em e 48c ( r ™) ffsc ) 
where E em « 2.8 • kTBB(? em ) is the average energy of a photon emitted by the blob, and 
O c (r ern ) = kT c (r em )/(m e c 2 ). The time t esc at which a photon emitted at r em escapes from the 



corona, may be estimated by t esc ~ t em + {r ou t — r em ) y N sc + 1 /c. The general trend resulting 
from this estimate is consistent with the observed approximate dependence At^i oc \n(E^/Ei) 
( Miyamoto et al. 1988 , Nowak et al. 199"8| ) of the hard time lags At^i on the mean photon energies 



Ek i in the energy channels k and /, respectively. 
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3. Numerical results and discussion 

We have performed a series of simulations with coronal parameters fixed to the ADAF 
self-similar profiles, but different assumptions on the blob temperature and its evolution. We find 
that the best agreement with observed features in the power spectrum and hard time lags of X-ray 
binaries is achieved if the blob heats according to Tbb cx r~ 1//4 . The power spectra of the X-ray 
signals in the different energy bands and hard time lags between adjacent energy bins resulting 
from a simulation with kTss = 0.2 keV at the last stable orbit of the accretion flow around a 
Schwarzschild black hole are shown in Figs. 1 and 2. 

The power spectra are characterized by a flat plateau below some break frequency and turn 
into a power law with index ~ 2 at higher frequencies. The break frequency increases with 
increasing photon energy, from <C 0.1 Hz at ^ 1 keV to ~ 2 Hz in the 25 - 100 keV band. This 
is in excellent agreement with the observed power spectrum of Cyg X-l flCui et al. 1997a| , 1997b, 



Nowak et al. 1998) which turns into a power law of index ~ 2 above a break frequency of a few 
Hz. No significant dependence of the high-frequency slope on the photon energy was observed by 
Nowak (1998) for Cyg X-l and by Motch et al. (1983) and Maejima et al. (1984) for GX 339-4. A 



flattening of the PSD with increasing photon energy, which is occasionally detected ( Brinkman et 



al. 1974 , Miyamoto and Kitamoto 1989j ) could be caused by the increase of the break frequency 



in the PSD, evident in Fig. 1. We point out that the simplistic scenario simulated here is only 
relevant to the high-frequency portion of the power spectrum, corresponding to variations on 
short time scales. On longer time scales (lower frequencies) the power spectrum in reality is 
dominated by the time series of multiple blobs which might correspond to accretion instabilities 
in the outer accretion disk. The observed broadband PSD is simply the product of the single-blob 
PSD calculated here and the power spectrum of the time series of blob recurrences. 

The hard time lags between adjacent energy bands are generally linearly proportional to the 
Fourier period, turning into a flat curve beyond a break period, which decreases with increasing 
photon energy. For the highest energy channels, the time lags depend on the Fourier period to a 
power less than 1. The maximum time lags between medium-energy X-ray bands (in the energy 
range accessible to the PCA onboard RXTE) are of the order ~ 0.1 s. This is in very good 
agreement with the hard time lags observed in Cyg X-l ( [Miyamoto et al. 1988 , Cui et al. 1997b 



Nowak et al. 1998| ). Even the flattening of the time lag curves towards At oc t 8 as observed by 



Crary et al. (1997) at high photon energies is well reproduced by our simulations. 

The photon spectrum produced by this scenario consists of a thermal hump peaking around 
~ 0.2 - 0.3 keV and a power law of photon index a ~ 1.5, typical of GBHC low-hard state 
spectra. Since we may assume that the global parameters of the corona and of the evolution 
and motion of the blob do not dramatically change on short time scales, the coherence function 



( Vaughan fc Nowak 1997 ) is close to 1 for most frequencies, except at very high frequencies, where 
photon counting noise will destroy the coherence, and at very low frequencies, where the stochastic 
sequence of blob excitation events will produce an incoherent variability pattern. Such a high 
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degree of coherence over a broad Fourier frequency range has indeed been observed for Cyg X-l 
(|Cui et al. 1997R |Nowak et al. 1998Q. 



When we increase the heating rate of the inward-moving blob we find basically the same 
trend as described above for the power spectra, but the hard time lags turn into much steeper 
curves, approaching At oc Pp (where Pp is the Fourier period) for a radial blob temperature 
dependence Tbb oc r~ 3 / 4 ( JShakura fc Sunyaev 1973 ). The temperature of the blob at the last 



stable orbit is tightly constrained by the observed soft excess in GBHC X-ray spectra, indicating a 
typical blackbody temperature around kTpB ~ 0.2 keV. If TsBi^in) differs significantly from this 
value, the resulting X-ray spectrum becomes inconsistent with the observed spectra of Cyg X-l 
and other GBHCs in the low — hard state. 

The maximum time lags and the turnover frequencies in the PSD and the time lag curves 
are determined by the drift time scale of the blob through the innermost, hottest part of the 
corona. All observables resulting from our model are insensitive to the value of the outer radius 
of the ADAF-like corona because the outer regions of the corona are inefficient in terms of 
Comptonization. With increasing drift velocity the maximum time lags become smaller, the 
turnover in the PSDs is shifted towards higher frequencies, and the turnover in the time lag curves 
toward a constant value occurs at shorter Fourier periods. If the radial drift time scale of the 
blob approaches the photon escape time scale for photons emitted near the black hole horizon, 
the resulting light curves become almost symmetric, and the hard time lags become small or even 
negative. Thus the choice of a radial drift velocity (3 r 1 %, much smaller than the free-fall 
velocity, is essential for our model to reproduce the observed trends in PSDs and hard time lags. 

A detailed parameter study, combined with detailed fitting of the observed X-ray spectra, 
PSDs, time lags and coherence functions of Cyg X-l and GX 339-4 will be presented in a 
forthcoming paper. 

In this Letter, we have only presented results from a single blob. However, in reality the 
observed X-ray signal will consist of the contributions from multiple blobs formed in an irregular 
sequence. Our results may be regarded as a transfer function tr^{t) describing the X-ray response 
of the corona in photon energy channel to a single blob moving inward. The effects of the 
convolution of a time sequence of elementary events with such a response function was briefly 
mentioned in Bottcher and Liang (1998) and demonstrated in detail in Nowak et al. (1999). Let 
b(t) be the number of identical blobs entering the corona per unit time. Then the observed signal 
in energy channel k will be fk(t) = (b * trk)(t), where * denotes the convolution. The Fourier 
transform of fk(t) is F^vp) = B(vp) Tr^^p), where capital letters denote the Fourier transforms 
of the respective functions. Since b(t) is independent of the photon energy, the resulting time lags 
are exactly given by the single-blob time lags as plotted in Fig. 2. The resulting power spectrum 
is PSDk(vp) = \B(up)\ 2 |Trfc(^ir)| 2 . The observed PSDs of GBHCs are typically given by a flat 
shape below a break frequency v\ ~ 0.1 Hz, a power-law PSD oc Vp for u\ <^ up <J V2 ~ 2 Hz, 



and a steeper power-law oc v a , often consistent with a ~ 2 for up ^ U2 (e. g., Cui et al. 1997c 



- 6 - 



Rutledge et al. 1998, Nowak et al. 1995, Grove et al. 1998| ). Thus, comparing to Fig. 1, we can 



deduce that the intrinsic power spectrum of the blob formation events, B(isp), must be roughly 
flat below v\, a power-law oc Up 1 for v\ <^ up z/ 2 , and again flat for vp ^ z/ 2 . This indicates 
that blob formation events happen on typical repetition time scales ~ 1 s ^ Ai rep ^ 10 s with an 
intrinsic power spectrum \B(uf)\ 2 oc Up. 

Another important point is the question of stability of the cool blobs since our simulations 
assume that the blobs are heating up slowly, but are maintaining their integrity. In the context 
of the inhomogeneous accretion equilibrium solution of Krolik (1998), this is equivalent to the 
condition that thermal conduction from the hot to the cold phase is negligible, which yields a 

— 1 3/2 

constraint on the minimum extent Rb of the cool blobs, namely Rb ^ 1.73 hr c Q c , where h 
is the scale height of the hot phase, r c is its half Thomson depth, and G c = ksT c / (m e (P) is its 
dimensionless temperature. However, in the case considered here, we only need to ensure the 
weaker condition that the time scale of conductive heating and thus cloud evaporation be longer 
than the drift time scale of the clouds through the innermost, hot region of the corona, which 
is ~ 1 s. A detailed time-dependent analysis of the stability of spherical clouds to evaporation, 
including the effects of viscous stresses, but neglecting radiative heating and cooling, has been 
presented by Draine & Giuliani (1984). If viscous stresses are unimportant, the fractional mass 
loss rate due to thermal conduction is 

M B _ 12 Kdass „ in _ 3 ( k B T c \ 5/2 o-2 -1 -i m 
M B --25k B Rln B ~ S ' W llOOkevJ ^ "» S (1) 

( Draine &: Giuliani 1984| ) where Mb and Mb are the total mass and mass loss rate of the blob, 



Kdass is the classical thermal conductivity ( Spitzer 1962]) , Rb = 10 7 i?7 cm is the blob radius, and 



tib = 10 20 ri2o cm 3 is the particle density in the blobs. Eq. (1) indicates that for ri2o Rj ^ 0.1 the 
blobs are expected to drift through the corona without significant mass-loss. 



4. Summary and conclusions 

We propose a new model for the rapid aperiodic variability of the X-ray emission of GBHCs, 
involving a cool, dense blob drifting inward through a hot, inhomogeneous Comptonizing corona 
and disappearing at the event horizon of a solar-mass black hole. We have demonstrated that the 
resulting spectral and temporal characteristics predicted by this model are in very good agreement 
with the observed X-ray spectra, power spectra, hard time lags and coherence functions of Cyg X-l 
and GX 339-4 if the hot corona has a radial temperature and density profile corresponding to an 
ADAF and the blob heats up slowly (Tbb oc r -1 / 4 ) as it moves inward. 

We have specified the conditions under which such cool blobs are stable against evaporation 
on the radial drift time scale of ~ 1 s, which is the case for 7120 R7 ^ 0.1. From a comparison 
with observed power spectra of GBHCs we can deduce that the power spectrum of blob formation 
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events must be given by a v F l power law, and that these events occur on a typical repetition time 
scale ~ 1 s <; At rep <J 10 s. 

Our model is in accord with the increasing evidence, from observational as well as theoretical 
results, that accretion onto Galactic black holes does not occur in a quasi-static equilibrium state, 
but that accretion instabilities lead to a rather erratic, inhomogeneous structure of inner accretion 
flows. The inward-motion of localized, efficiently radiating cool blobs through an inhomogeneous, 
hot corona then provides a very natural explanation for the observed hard time lags in the 
aperiodic variability of GBHCs. 

We thank the anonymous referee vor very helpful and constructive comments. This work is 
partially supported by NASA grant NAG 5-4055. 
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Fig. 1. — Power spectra of light curves in five photon energy bins for a blob drifting inward 
with P r = 5 ■ 10 -3 through an ADAF-like corona and heating up according to Tbb °c r -1 / 4 , 



kTBB(ri n ) = 0.2 keV. Coronal parameters: kT c (ri n ) 
r in = 10 7 cm, r out = 10 9 cm. 
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Fig. 2. — Hard time lags between adjacent photon energy bins as a function of Fourier period Pp; 
for parameters see caption of Fig. 1. Thin curves are the theoretical continuation of the phase lag 
curves at long Fourier periods which are not covered by our simulations. 



